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Abstract 

We study the Schwinger model on a lattice consisting of zeros of the Hermite polynomials 
that incorporates a lattice derivative and a discrete Fourier transform with many proper- 
ties. Such a lattice produces a Klein-Gordon equation for the boson field and the exact 
value of the mass in the asymptotic limit if the boundaries are not taken into account. 
On the contrary, if the lattice is considered with boundaries new modes appear due to 
aliasing effects. In the continuum limit, however, this lattice yields also a Klein-Gordon 
equation with a reduced mass. 
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1 Introduction 



Since lattice theory is intended as a nonperturbative way to solve problems in Quantum 
Field Theory both from analytical and numerical points of view, it is important to study 
finite lattice models in order to support consistent numerical calculations as well as the 
belief that an infinite lattice is in fact the continuum limit of a finite one. The Schwinger 
model has proved to be an excellent theoretical laboratory to study physical properties 
of fermions on the lattice [TJ |2l El IH IS]- Several approaches and derivatives has been 
used to study the mass spectrum, chirality, charge screening and the chiral condensate, 
among others. In fact, in jl] it has shown that different types of fermions derivatives, 
namely the naive, the Wilson's and the so called modified SLAC derivative work well for 
the Schwinger model, yielding the correct boson mass in the continuum limit, but also 
producing a infinitely heavy boson mode. In the past years, two different differentiation 
matrices (nonlocal lattice derivatives) have been used to circumvent the restrictions of 
the no-go Nielsen-Ninomiya theorem [HI Ej . One of them is the Trigonometric derivative 
defined for 27r-periodic functions 0. The other, which can be called Hermite derivative, 
is suitable for nonperiodic functions decreasing rapidly to zero at infinity. Each derivative 
has a well-defined discrete Fourier transform. If the number of lattice points is an odd 
integer, the spectrum of the Trigonometric derivative coincides with that obtained by 
using the SLAC derivative. 

In this paper we study the lattice Schwinger model using a lattice consisting of zeros 
of the Hermite polynomials that incorporates a lattice derivative and a discrete Fourier 
transform. We obtain the correct boson mass in the asymptotic limit if the boundaries 
are not taken into account, however, new boson modes appear in the case in which the 
boundaries of the lattice are considered. This effect is completely due to the aliasing 
inherent in a finite lattice. 

Let us start by considering the explicit form of the derivative used in this paper. The 
N X N Hermite differentiation matrix D is given hj D = SDS~^ where 



S" is a diagonal matrix ^ which nonzero elements are given by exp[—xl/2]H'^{xk) , k = 
1,2, . . . , N, and Xk is the kth zero of the Hermite polynomial Hn{x). The prime means 
differentiation with respect to the argument. This matrix yields exact derivatives at the 
lattice points x^ for functions of the form exp(— x^/2)/(x), where f{x) is a polynomial 
of degree less than A^, and in general gives approximated values to the derivative of a 
differentiable square-integrable function, since a function like this can be expanded in a 
Hermite series. Taking into account that 
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the elements of the matrices D and D coincides asymptotically up to a factor (— l)-^"*"*^ 
and therefore, in this limit D becomes a skew-symmetric matrix with a simple structure. 
The differentiation matrix is diagonalized by the discrete Fourier transform F. Thus, we 
have an equation of the form F^DF — iP where P is a diagonal matrix containing the 
eigenvalues of —iD. The discrete Fourier transform is given by 



N-l 



^Ok = Pkk' = Pkhk', 

1=0 

where 



Here, Xk and pk are roots of HnH), and if k is fixed, they have the asymptotic form 

Xk^Pk^ -^j^ k, k integer. (2) 

In this hmit, the discrete Fourier transform behaves as Fj^ — (AxAp/27r) 6*^^^*= + 
0{1/N) where Ax = Xk+i — Xk — Ap = pk+i — Pk — 7r/y/2N. Since both Aa; and Ap are 
dimensionless and equal, we can write this asymptotic expression in terms of the Riemann 
measure A^ = Ax = Ap 

Fjk = ^ e^^^P* + 0{1/N). (3) 



/27r 

2 The lattice model 

Our study of the lattice Schwinger model is based on the Heisenberg equations and canon- 
ical quantization, therefore, lattice artifacts will be used only for the space variable since 
time can be considered to remain continuous. In the temporal gauge Aq — 0, the Hamil- 
tonian of the Schwinger model is given by 



H ^ dx 



^E'^{x) + iip^'a^ {dx + ieA{x)) ip{x) 



(4) 



where A{x) stands for the component A\{x) of the gauge field. To get a discretized version 
of this Hamiltonian it is necessary to consider the action of any differentiation matrix D 
on the fermionic variable in more detail. Since 



N 



{Dr)i^y^p>^^rk. 



k=l 



where "0" is the vector whose elements are the values of the a-component of the fermionic 
field iIj{x) at the lattice points and time t, the discretized kinetic terms become nonlocal 
for a finite number of nodes. Thus, for finite N the covariant derivative must be replaced 
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by a matrix differentiation with links. This task can be carried out by noting that under 
the transformation Xk —>■ Xj the fermionic field transforms as 

^{xk) ^ e"^-7 ^'^^^''^^xj) = e'^'^e-'^'^^^ixj), (5) 

where we have written J^'' A{^)dC, = bk—bj. Therefore, the covariant derivative d^+ieA^x) 
is discretized in our scheme as 

V = t/t o D o f/, (6) 

where U is the unitary matrix with elements Ujk = 5jkG~^'^^^ and o denotes the element- 
wise matrix product. Since the differentiation matrix (^J) become local derivative in the 
continuum limit, V yields the correct form of the covariant derivative. 
Thus, our lattice version of (j^ is given by 

HL = ]^E'E + iij\a-,®S/)ij, (7) 

where E and are the vectors constructed with the values of the fields E{x) and ip{x) at 
the lattice points and the superscript t means transpose. 

Taking into account the similarity transformation (jHl), Hl can be written in the form of 
a free Hamiltonian 

HL = Hp + H^ = ]^p'D-''p + i¥{a,®D)'^, (8) 

where = Uip"' and the discretized form of Gauss' law 

p = DE = e^'^^ o = e^'^^ o 

has been used. In order to have a nonsingular matrix D, the number of nodes N must be 
an even integer. The rth value of the axial current j is given by 

and the fermion fields tp, tp"^ satisfy the equal-time anticommutation relations 

{i;l'^,^^} = 5rs5ap. (9) 

Note that in order to make the connection to the continuum case, should be multiplied 
by to obtain the Riemann sum of (0]), whereas should be divided by to obtain 
a delta function. 

3 An equation for the charge density 

To find an equation for the bosonic variable p on the lattice we begin by considering the 
lattice version of dop = [p,H]/i, which can be written as {dop)r = [Pr,HL]/i. The error 
in the derivative of a field produced by the use of a differentiation matrix is negligible if 
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the number of points is great, therefore, we assume a lattice with a sufficiently great 
number of nodes. There are, however, two possibilities to consider in this asymptotic 
limit. One is what can be called an infinite lattice, i.e., a lattice with no boundaries and 
therefore, no border effects. The other is a finite lattice with definite boundaries. The 
following first calculations are common for both cases since they only use the fact that 
is sufficiently great. 

By using Q and the skew-symmetry of D we obtain the continuity equation 

= e (^tv,^(D^/^), + (D^t")^^^^^^) (10) 

Now, let us compute 

{dlp\ = -[iDj)r.,HL]/t = -[{Dj\.,Hp]/t - [{D3\,H^]/i. (11) 
Proceeding as before, we find 

[{Dj\, H^]/i = -e{^l[{l2 ® D^)%r + 2[(l2 ® D)m[{l2 ® D)%, + [(I2 ® D^)^^,) 

= -(^V)r, (12) 

where I2 is the identity matrix of dimension 2. Therefore, ()11|) becomes 

{dlp)r - {D'p)r = -^[mr,p'D-'p]. (13) 

To obtain the asymptotic form of the mass term it is necessary to compute the Schwinger 
commutator [jj,Pj']. To this end we compute the commutator of Fourier transforms 
[jk,Pk']- First, we expand in its Fourier components 

k 

where F stands for the discrete Fourier transform (j21). Thus, the Fourier representation 
of pg = e^t^g is 

PQ = eJ2\-^i^k'FlgFgk' + bkblFgkFlg], 

k,k' 

where we have used the orthogonal property of the spinors. Therefore, the Fourier trans- 
form pk = Y,q FkqPq bccomes 

Pk = eY. [aW'Fl,^F,,.Fl + h'.hl,F,k„Fl^Fl]. (14) 

q,k',k" 

Note that the sum of products of elements of the discrete Fourier transform F involved 
in this equation are proportional to the generic asymptotic form 



i-5^e-^^Ae, (15) 



27r 
1 



where P stands for a sum of three lattice momenta, for instance, P = pk' — Pk — Pk"- It 
is worth to notice that a sum hke (jl5p is exactly a Kronecker delta if P is the sum of two 
lattice momenta (since F is unitary). Therefore, (jl5|) becomes an exact delta whenever the 
sum of two of the three momenta summing up P is again an allowed lattice momentum. 
But this is not always the case, even for uniformly spaced momenta. However, we can 
consider an infinite lattice, i.e., a lattice with a very great number of nodes with no 
boundaries and extending to infinity. Since there is a momentum for each lattice site, the 
lattice momenta do not have neither a maximum nor a minimum and therefore, the sum 
of two evenly spaced momenta is again an allowed momentum, yielding that 

^^e-,PAe-<5(P)Ae (16) 
q 

In this case, (fT^ becomes 

Pk = e^= ^[alttk'+k + bk'bl,_k]. (17) 

A similar computation can be done for jk- To obtain the value of [jk, Pk'], it is necessary to 
write this commutator in terms of operators that are normal-ordered with respect to the 
vacuum state filled with particles of any negative energy and antiparticles of any positive 
energy, defined as 

ivac)=n4n^iio)- 

fc<0 fc>0 

Thus, after the use of the anticommutation relations 

{al,ak'} = 6kk', {bl,bk'} = 6kk', (18) 
the normal-ordered commutator becomes^ 

[jk,Pk'] = y^(QmQ-m+fc+fc' - oL-fc'^m+fc' + b^m-k^k'^rn - b^rn-k'^rn+k) ■, 

m 

and, as it can be seen, it annihilates the vacuum state unless k' = —k. Therefore, 

A,^ 

[jk,Pk'] = e^6k',-k^ '^{(^L^^m - a^^^^am-k + ^L^m - &L+fc^m+fc) 

m 

and the action of this commutator on |vac) is given by 

[jfc,Pfe']|vac) = A^ A;5fc/ _fc|vac). 

TT 



""^ Again, in order to make the connection to the continuum case, H18|) should be divided by to 
obtain a delta function. 
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Thus, the inverse transform of [jk,Pk'], is 

2 2 

[JvPA = - V E kF,kFr,-k = - ^ E '"F^'^^lf (19) 

TT — ' TT ^ — ' 



The product kFjk is the asymptotic form of pkFjk = — i J2i DjiFik, as it can be seen 
from (0)- By using the fact that FF^ = 1, we obtain finally the Schwinger commutator 
on an infinite lattice 



.e2 



[jj,Pr]=i—Djj,. (20) 



vr 



Therefore, the right-hand side of (fT^ becomes the correct value of the mass term 

m^Pr = hiDj)r,p'D-''p] = -pr. (21) 
Zt 71 

Let us consider now the case of a finite lattice, i.e., a lattice with a great number of 
nodes and definite boundaries. Since momenta and sites have same numerical values, the 
allowed momenta lie in the same interval as the sites do. Thus, in the case of a finite 
lattice, the sum of evenly spaced lattice momenta may lie outside the interval of momenta 
and the exponential in (jl5|) produces aliasing in the momenta domain^. To explain in 
more detail the aliasing effect let us use (0) to rewrite (fTH|l as 

^ E^^P(^^^^)^^ = ^ E^^P[^27rg(7r/s:/4)]Ae, (22) 

9 9 

where k, = K/N, g is a nonzero integer of [—N/2, N/2] and K is a sum of integers of 
[—N/2, N/2]. In the asymptotic limit k becomes dense in the real line and two momenta 
Pk and pk' are indistinguishable for the model if their indexes, k and k', are congruent 
module 4/7r, i.e., if they yield the same sum ()22|). In other words, a momentum pk is 
distinguishable if its index k lies inbound in the interval [— 2/7r, 2/7r]. Note that this is the 
interval of periodicity of the exponential function of ()22|) . Thus, the sum of two lattice 
momenta (whose indexes are in [—1/2, 1/2]) can be a lattice momentum again with k lying 
in [—1/2,1/2], or a totally new mode with k, G [— 2/7r, — 1/2] U [l/2,2/7r], or an aliased 
and repeated mode with k e [-1, 1/2 - A/n] U [A/n - 1/2, 1]. The interval [4/tt - 1/2, 1] 
is aliased with [—1/2, 1 — A/n] and [—1, 1/2 — A/n] with [A/n — 1, 1/2], such as depicted 
in Figure 1. 



new 



new aliased 




repeated 



t -2/71 -7/2^ 
1/2-4/n 



lattice momenta 



m 2/n t 

4/71-7/2 



7 k 



^Aliasing is an effect associated with discretization of periodic continuous systems. In a video, for 
instance, the wheels of a car appear to be static or to rotate in reverse. This effect is due to ahasing of 



the wheel rotation frequency (see for example P|) 



7 



Figure 1: Aliased modes appearing in a finite lattice. 
Since N is sufficiently great, (fTH|) still holds in the form 



where pk + Pk" stands for the value (aliased or not) corresponding to the sum pk +Pk"- 
Note that the aliasing does not occur if one of the two momenta has an index \k\ < 
4/iT — 1 ^ 0.27, i.e, k G [1 — 4/7r,4/7r — 1]. For other values of k, the aliasing effect 
changes the charge density (jl7|) in the momenta space to 

Pk' = e^= aJa^+K' + ^ aja^+^,_i 



whereas the axial current becomes 



As before, the normal-ordered Schwinger commutator [j^jPfc'] annihilates the vacuum 
state unless k' = —k (or /c' = — /c). In this case the action on |vac) is 



- k |vac), < ^ — 1, 

[jfc,p,,]|vac) = -Ae5fc,_fc<;iV(l-i)|vac), < 1 - ^, (24) 

_iV(i -l)|vac), «:>^-l, 



e2 



TT 



and now (|19|) becomes 



TT — ' 

|A:|<(4/7r-l)Ar 

+ z— A^2(4/^_ J- sinpfe(xj- -Xj-/)A^. 

A;>(4/7r-l)Af 

The ffist sum of the right-hand side corresponds to the values of the momentum index 
for which aliasing does not occur and it approaches to (f^n|) in the asjTiiptotic limit. The 
second sum can be approximated by an integral to yield 

2 2 

[j„ p,,] = D,y + lU^I-K - l)(-l)^+^'5(i/,,0 (25) 
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where S{y) = [cos(4/7r — l)y — cosy/2], i/jj' = 'k^N 12 {xj — Xji) and we have used the 
fact that — Xk) = {—l)^^''Djk. The substitution of this commutator in the mass 

term of (jl3p yields two terms, one of them corresponds to the first term of the right-hand 
side of (123) and it gives the correct mass e^/vr. Due to S{yjj'), the remaining term yields 
a mixing of modes pr' at locations different from r for a finite [cf. Eq. ()21|)]. However, 
since \S{yjj')\ < 2, the behavior of S{yjf) Djj> is dominated by Djj/ and in the continuum 
limit, i.e., in the limit iV — > cx), S{yjj') Djj/ becomes essentially Djj/ due to the singular 
behavior of the matrix differentiation [7j. Thus, in such a limit, the main contribution 
of S{yjj/) is given by its first maximum (see Figure 2) located at yjj/ ~ ±5.14, i.e., at 
Xj — Xj/ ~ ±7r/-\/2iV = ±A^, or equivalently, at j' = j ±1. 

m 

10 20 30 40 50 

y 

Figure 2: Plot of the function S'^{y)/y'^ 

Since |S'(5.14)| ^ 1, we have that \S{iijj/) \ Djj/ ^ Djj/, N ^ oo, and therefore, in the 
continuum limit (^^j) behaves as 

2 2 

[j„ p,/] = D,y - ^^(4/7r - 1) D,y, (26) 

yielding 

= (^ - -)e^ = O.lSe^ (27) 

2 TT 

for the square of the mass of the Schwinger model on a lattice with boundaries. 
It is important to stress that aliasing is not just a feature of the discrete Fourier transform 
used in this paper and that aliasing effects should appear in lattice field theory whenever 
a finite lattice and the usual discrete Fourier transform are used to compute a quantity 
which depends on the sum of more than two momenta, as in Eq (jSHI)- The equivalent 
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form of this equation (using the usual Fourier transform) is 



l2Tl-^-{k' -k-k") 

e \ 

where fc, k' , k" = 1,2, . . . , N, and we have omitted the normahzation factors. In this case, 
we have again ahasing effects since the sum k + k" may be an integer greater than N and 
therefore may he outside the finite lattice. Quantities involving more than two Fourier 
transforms appear for example, when computing vertex functions. 



4 Concluding remarks 

We have introduced a nonperiodic lattice that incorporates a lattice derivative and a 
discrete Fourier transform with useful properties. Such a lattice has potential use in the 
study of quantum-field problems. In the case of the Schwinger model, the use of this 
nonperiodic lattice yields the correct bosonic mass if the lattice is considered unbounded. 
On the other hand, if the lattice has boundaries, a great number of new modes appears 
due to the inherent aliasing in the finite Fourier transform. In fact we have shown that 
in this case, the boson acquires a smaller mass. Since aliasing is an effect present in the 
usual discrete Fourier transform, the appearance of new modes or aliasing effects should 
be expected in any lattice approach where the standard discrete Fourier transform is 
used. Such effects, mainly due to the boundedness of the lattice, indicate that the results 
obtained by numerical calculations using (necessarily) finite lattice models of quantum- 
field problems defined on unbounded manifolds may not be completely correct if the 
aliasing effects are not taken into account. 
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